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ABSTRACT 

A low energy theorem concerning a matrix element of the QCD axial anomaly is 
tested against different instanton models. In the chiral limit the theorem is ful¬ 
filled by the Diakonov&Petrov model, whereas it is violated by single-instanton 
approximations. Beyond the chiral limit the theorem, and two relations estab¬ 
lished for other matrix elements of the QCD axial anomaly, result to be violated 
also by the Diakonov&Petrov model. 


1. Introduction 


In quantum field theory the symmetry of the classical lagrangian may be destroyed 
by quantum correction. In gauge theories the most important examples of this kind 
are the axial anomalies in Electro Weak theory and in QCD. The axial anomaly arises 
from noninvariance of the fermionic measure against axi^ transformations in the path 
integrals of the theoryhl (see also higher-loop correctionsH) . In Euclidean QCD-I-QED 
the axial anomaly reads 

aV; = ^ Q]FF + 2 ^ (1) 


where = S/V’/T/iTst/’/(/ = u,d,s) is the quark singlet axial current, -0/ the 
quark field, Nf the number of the flavors, g the QCD coupling constant, 2GG = 
Qa^ gluou field Strength operator, Nc the number of colours, e the 
QED coupling constant, Qf the electric fractional charges of the quarks, the 
photon field strength operator. We have explicitly included the contributions of the 
current masses of light quarks ruf. It was pointed outB that this equation leads to a 
nontrivial low-energy theorem {LET) and to useful relations. 

In the present paper we consider the matrix element of eq. O {i) between vacuum 
and two-ph<^ons states and {ii) between vacuum and meson states. 

Analysis^ of case (i), and nonvanishing of the r] meson mass even in chiral limit 
due to axial anomaly, leads to the conclusion that in the kinematical configuration of 
zero virtualities of each line in the form factor the matrix element of the divergence 
of the singlet axial current vanishes, implying the following LET : 


<0A/^GG|27> = +2i^m;(0|q75V>/|27), (2) 
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where — ei^uqi,fi and qi, e* {i = 2,3) are respectively the momenta and 

polarization vectors of photons. 

Eq. (^) gives the exact answer for the matrix element of the gluonic operator. Of 
conrse gluons can couple to photons only by virtual quarks. In perturbation theory 
the left side of the eq. (|^) is at least O^g'^e^), whereas the right side of this equation 
does not contain any strong coupling at all. It is evident that the LET can be 
fulhlled only out of the framework of perturbation theory. The factor at the left 
side must be completely cancelled by som^nonperturbative contribution, such as the 
one, of order g~'^, provided by instantonsB. The computation of the left side of eq. 
(1^) amounts to calculating the instanton contribution in the Euclidean space to the 
three-point function 

r^u{xi,X2,X3) = {0\T[g^GG{xi)f^"'{x2)jt"'{x3)]\0), (3) 

UfT' being the electromagnetic current) which contains needed information on the 

matrix element ( 0 |GG| 27 ) and also useful physical information concerning r/ —>■ 27 
decay, 27 - transitions in heavy quarkonia etc. 

The matrix elements of eq. O between vacuum and meson states lead to the 
relations 

(0|iV/Y|^GG|??) = 2zm,(0|V’l75^,|7), (4) 

^ - mrf)(O|'0V375^/>|7r°), (5) 

which provide, together with LET, more stringent tests for instanton models. 

Instantons - whose presence in QCD is a well established fact, nt least at a phe¬ 
nomenological level and in numerical simulations of QCD vacuumO - constitute the 
main input to our calculations, especially in connection with quark propagation. To 
our present knowledge the instanton structure of QCD vacuum is coiicentrated in an 
average size p and in an average interinstanton distance R, such thatOQ 

p= 1/3/m, R = lfm. (6) 

Therefore the packing parameter [p/R^ = 0.012 is small, legitimating independent 
averaging over positions and orientations of the instantons; moreover in sprue cases 
it is even possible to usp^the quark propagator in the single-instanton £eldlI3. 

In ^mrevious workEdI LET was used as a test for the Diakonov-Petrov {DP) 
ansatzlldiy of the low-energy QCD effective action in the chiral limit. The anzatz, 
based on an interpolation formula for the quark propagator in the held of a single 
instanton, does satisfy LET to ~ 17% accuracy. 

One of the aims of our present work is to formulate more stringent tests of the DP 
model, by taking explicitly into account the contribution of the current quark masses 
ruf. We also test the validity of different instanton vacuum models by calculating 
in the range of small virtualities and testing it against LET. 

In section 2 we calculate in single-instanton approximation(S'/%). Irnsection 
3 we rederive the DP effective action starting from results of Lee&Bardeenlldl (LB) 
and accounting for current quark masses mj. In section 4 we calculate by this 




effective action and check it by LET. In section 5 we calculate in the same framework 
both sides of relations and and compare each left side with its respective right 
side. 

2. Single-instanton approximation 

A lot of papers have been devoted to instanton calculus to single-instanton ap- 
pmmoation (see for example the calculations of the correlator of two vector currents 
inllSUZiy). In this approximation the calculation of amounts to computing integrals 
like 

J d^z±dU±{T{g^GG{x,)f;^{x2)jr{xs))± = 

[ d‘^z±{GG{xi)±Tr{'y^S±{x2, X3)'yuS±{x3, X 2 )), 
f 

where z± and U± are the position and orientation of the (anti-)instanton, assuming 
instanton integration sizes to be peaked at p, eq. (^). Here 

(/GG(x))± = ±/(x - j) = ±—(8) 

[p^ (x - zy\ 

and S± = {iD± + im)~^ is the full quark propagator in presence of a single 
(anti)instanton. 

Starting from the expression of S± given in ref.ty, we obtain 


T^,u{Xl,X2,X3) = 


N 192iVcp® 


X 


V 

.h(x 2 - z) + h(x 3 - z) 

(X2 -X3)2 


J d^zE^{xi - z) 


h{x 2 - z)h{x 3 - z) 
{X2 - XsY 


h{x 2 - z)h{x 3 - z)]e^,yap {X 2 - z)o,{x 3 - z)p, 


where h[x) = (p^ -|- x^) Then the correlator, that is the Fourier transform of 
the three-point function to be tested against LET, results to be 




( 10 ) 


rl rl 

x4! / ada db{a{l - a + ab{l - -x- 

Jo Jo OP2,aOP3,P 

as can be shown by applying the Feynman integration technique. We have set 

fiql) = J d^xi exp(*gixi), 

J{P) = J d^Y [Y^ + expiPy, 


( 11 ) 









having introduced the 8-dim vectors P = (^2,^3), ^ = (2/2,1/3), with 

q2 + bq3 gs 2 2n M 

The latter integral (0) can be easily calculated and is reduced to the MacDonald 
function, i. e., 

J(P) = ~K,(Pr), P = (Py/\ (12) 

In the limit of small momenta 

1 Pr Pr 

K,{Pr) = + + 1/2] + (13) 

where C = 0.577215... is the Euler constant. 

Eqs. (pUP, (]T^, ( ]T3|) imply that T^u^qi, q 2 , qs) is divergent for qf 0, that is, the 
model badly violates LET. Nor does agree with LET the improved single-instanton 
approximations, where the zero-mode contribution to the propagator has been mod- 
ihed by replacing m —> m*, where m* is the effective mass of the quark accounting 
for the influence of the surrounding instantons). We conclude that single-instanton 
approximation badly violates LET, therefore it is not suitable for calculating three- 
point functions at small momenta. By the way we note that two-point functions of 
vector currents are not so sensitive to large distance effects. The violation of LET 
- and in particular the divergence in massless quark aproximation - is related to the 
slow decrease of at large distances. Indeed we have neglected rescattering effects 
of quarks by other instantons, which, during quark propagation, lead to the formation 
of the constituent quark, producing a suitable effective mass and providing needed 
exponential decrease with distance. Such effects can be described by an effective 
action. 

3. A low energy QCD effective action 

It is natural to choose the singular gauge for the instantons in describing many 
instanton effects in the propagation of the quarks. In the case of a small packing 
parameter it is possible to do the following ansatz for the background instanton held: 

N+ V_ 

H A) + Z! A), (^± = (^±, U±, P±)), (14) 

+ 

where, Zi, Ui and pi the position, orientation and size of the i-th instanton. The canon¬ 
ical partition function of the A+instantons and N_ antiinstantons can be schemati¬ 
cally written as 

Zn+,N- = / detArexp(-I4) J]] ,AzidUidn{pi), (15) 

i 

where Vg is the instanton-(anti)instanton interaction potential generated by the gluon 
held action and detj^ is a quark determinant in the instanton held. The main assump¬ 
tion of the instanton model is that Vg is repulsive at small distances between instanton 




and antiinstanton. This should provide the stabilization of the instanton sizes and of 
the interinstanton distances, as discussed in the introduction. We mainly deal with 
det^, which describes the influence of light quarks. 

Lee&Bardeenh3 (LB) calculated the quark propagator in a more sohsticated ap¬ 
proximation than SIA. In particular they found that 

detN = detB, Bij = imdij -|- (16) 

where aij is the overlapping matrix element of the quark zero-modes generated 
by instantons. This matrix element is nonzero only between instantons and antiin- 
stantons (and vice versa) due to the chiral factor in ) i- 6-) 

a_+ = - < <h_,o|ia|$+,o > . (17) 

The overlapping of the quark zero-modes causes quark jumping from one instanton 
to another one during propagation. 

It is clear from ([T6|) that for iV+ ^ iV_ detjq ~ so the fluctuations of 

\N+ — NB are strongly suppressed due to presence of light quarks. Therefore we 
assume N+ = = N/2. 

Following the procedure suggested in ref.h3, we get the fermionization of LS’s 
result, i.e., 

detisf = J exp{J d'^x 'i/jlid'i/jf) 

N+ nJ (18) 

X ,'ipf])), 

f + 

where 


Eq. (p!8D coincides with the ansatz for the hxed N partition function postulated by 
DP, except for the sign in front of V±. Keeping in mind the low density of the 
instanton media, which allows independent averaging over positions and orientations 
of the instantons, eq. (PIBf ) leads to the partition function 

Zn = J D'lpD'ip'^ exp{J d'^xipUd'ilj)W^^ W^~, ( 20 ) 

where 


W±= I d^^±Y[{V±['i/j}7pf] +imf) = 



and 

r d^kdP 

J±{z)fg = J -j^;^exp{-i{k-l)z)F{k‘^)F{f)ij}{k)l{l±'y 5 )'ipg{l). (22) 





The form factor F is related to the zero-mode wave function in momentum space 
aiid is equal to 


F(k'‘) = [/o(«)A'„(() - /i(()A'.(«)|, t = 


(23) 


4. Calculation of the correlator by the DP effective action 

In quasiclassical (saddle point) approximation any gluon operator derives its main 
contribution from instanton background. In the following the operator g^GG{x) will 
be considered. Owing to the low density of the instanton medium, it is possible 
to neglect the overlap of the fields of different instantons. In that case, the matrix 
element of g‘^GG{x) with any other quark operator Q is 


{g^GG{x)Q)N = J D'lpD'ip^ exp{J d'^x'ipHd'ip) 

X {n+ {Waa+(x)Q) W^*-^w!!- + N. {Waa.(x)Q) , 


(24) 


where 



X 1 / .r / X .m Nr. 

z)det(J±(s)+jj^). 


(25) 


It is useful to introduce the external fields k{x), coupled to g^GG, and a, such that 
D = d — zeQ/O— Starting from we find that the partition function Z[k, a] de¬ 

scribing mesonstll in presence of such external fields is 


Z[K,a] = J Di+D<i.exp(-Wl<t+,<i.]), 


(26) 


where 

hT[<I)+, <I)_] = J d'^xiWa + Wb- Wr), 

= (A'/ - l)^(n V7‘<ief4+)<'''-‘)-‘ + (+ ^ -), 


Wb = 

Wr = 

A± = 


Nr 




-Tr{m ($+ <!>-)), 


y^Trln 

f 


id + iF^ ($+24++ <I)_24_) 
id + iruf 




2(1 ± 75 )- 


( 27 ) 









The saddle point of the integral (p6D is located at {^±)fg = MfSjg, a self-consistency 
equation for the effective quark mass, i. e., 




d^k M]F\e) 

(27r)4 M|F4(F) + P 


= N + 


rrifMfVNc 


(28) 


being imposed. Eq. (^) describes also the shift of the effective mass of the quark Mf 
due to current mass mj. Expanding (^8]) with respect to mj, we have Mf = Mo+'ymf, 
where 


-1 


roo 

.2 / 112 


'' ' ''To {M§F*{k^) + 


(29) 


Such integrals converge due to the form factor F{k‘^). Assuming for the instanton 
model parameters p and R the values (Q) - which correspond to Mq = 340 MeV - we 
find 

7 = 2.4 . (30) 

The quark condensate is then given by 


^.-ly-idZN 


/^=0 


N^Mq ^ -{2QhMeVf r^N, ^ 


27r^p2 




(31) 


This quantity can be also calculated by formulaS —i < >Euciid= iTr S, yielding 

• , ..r f MoF\k^) 

—l < >Euclid — —4:Nc J 


{2Ey M^F^{k^) + k^ 


= -(255 MeE)^ ~ 


(32) 


pR^' 


Although coming from different formulas, predictions (^) and (^^ have the same 
parametric dependence and agree with the phenomenological value, i. e., 

-z < V'V >Euciid= -(240 - 250MeE)l 
The three-point function can be derived from the functional relation 


rgu(,Xi,X2,X3) = 


5W[k, a] 


K,,a=0’) 


where0 


W[K,a\ = Y^Trln 
f 


6K{xi)6ag,{x2)6a,y{x3) 

W + iMfF^ (A+ + A_) 


id + inif 


(33) 


(34) 


It is clear from eq. ([3^ ) that the vertex factors in the diagram of the process are 
eQfjg, and iMffF'^Nj^'y^. Taking the Fourier transform of (|33|) , we get 




(?i, ?2, gs) = f{ql)Nce^ f{ql, ql qj), 

f 


( 35 ) 













where Tf{q‘f,q 2 ,ql), the form factor coming from the diagram of the process consid¬ 
ered, may be calculated analytically if we approximate the form factor F by 1 . As a 
result, the left side of eq. (^) reduces to 

( 36 ) 

which coincides with the hrst term at the right side of eq. (|^). If we take into account 
the form factor F in ([35|), and give the model parameters the values (^, in the chiral 
limit we hndllU a variation of ~ 17%. Beyond the chiral limit the left side of (^) 
receives the contribution 


N 

- Y. (37) 

This quantity has to be compared with the explicit contribution of the current quark 
masses to the right side of (H), which, too, may be calculated by formulas (|3^ and 
(|33|) substituting by 2imf in the 75 vertex. Approximating again F ~ 1 , 

we obtain 

2*E"^/(0|V'/75^/’/|27) = (38) 

/ f ^ 

The ratio of ( 0 ) to (^) at the parameter values (^) results to be 

- 0.27pMo = -0.28 (39) 

and not 1, as demanded by LET. We stress the neat contradiction with the theorem, 
not only in magnitude but also in the sign. So the DP model (eq. (^)) fails to 
reproduce LET beyond chiral limit. 

5. Other refined tests for instanton models 

The matrix elements of the anomaly (|l|) between vacuum and 7 -meson or 7 r° states 
lead to more stringent tests of the DP model. 

The partiti^ function ( 0 ) describes mesons with matrices <I>±, whose usual de¬ 
composition is^ 

<l>± = exp(±^0)M(Texp(±^0), (40) 

0 and a being Nf x Nf matrices. The saddle-point condition demands cr = 1 , <I)± = 0. 
The usual decomposition for the pseudoscalar helds 0 = ^^7 be used, where 

TrXiXj = 26ij and 08 ( 3 ) can be identihed with the 7(7r°)-meson state. 

Firstly we consider the matrix element in which the 7 -meson is^volved. Ne¬ 
glecting the small admixture factor (~ 0 . 1 ) with the pure singlet stateEH), the matrix 
element of the divergence (H) between the 7 -meson and the vacuum leads to eq. (P, 
which can be used as a test for instanton models. 






As it is clear from previous considerations, the factor g^GG generates the vertex 
iM and the //-meson gives rise to iMXsF^'j^. The structure of the mass 
matrix is 

M = Mo + 7(m«(- - + rud —^). 


Then at small q 




d^k MoF\k‘^) 


(27r)4Mo2F4(F) + F' 
Applying the same procedure that led to eq. (p8[), we get 

~ 2 N 


The right side of (0) is 


2m^(0|//^j75V>s|7) = 2m^(-77)4A^< 


On the other hand eq. 


v/3^ 
yields 


d^k MoF\k^) 
(27r)4 M^F^k^) + P 


m.o 


pR^ ‘ 


2m,{0\i>lj!ii>,\p} = 2m,( - —)i < > 


(41) 


(42) 


(43) 


(44) 


Now let us confront such equations with some consequences of the relation (§), where 
the TT^-meson is involved. Repeating for this case the calculations applied to relation 
(^), the left side of equation (|) yields 




(45) 


while the right side results in 


2z(m„ - mrf)(O|//’'''y75'0|7r°) = 2(m„ - m^)/ < tjj^ > . 


The ratio of 
d) yields 


(46) 

to (^41) equals the ratio of (^5|) to (|46|) and at the parameter values 
N 

- = 0.66 (47) 

and not 1, as demanded by relations (^) and (|^). Again the DP model (eq. (^OD) 
strongly contradicts some consequences of the operator equation (|^ beyond the chiral 
limit. 


6. Conclusions 

We recall the main results of our treatment. LET and relations (^-(K^piistitute 
sensitive tests for instanton models. The single-instanton approximationt3,113jl2 and 

















its modificationsLl,B badly violate LET. On the contrary the DP model, which takes 
into acconnt the multiinstanton effects in the propagation of the qnarks, agrees sat¬ 
isfactory with this theorem in the chiral limit. However, beyond such a limit, even 
this model violates LET and the relations (§)-(§)• 
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